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A potential thermoelectric material CuA102 is theoretically studied. We first construct a model 
Hamiltonian of CuA102 based on the first principles band calculation, and calculate the Seebeck co- 
efficient. Then, we compare the model with that of a well-known thermoelectric material NaccCo02, 
and discuss the similarities and the differences. It is found that the two materials are similar from 
an electronic structure viewpoint in that they have a peculiar pudding-mold type band shape, which 
is advantageous as thermoelectric materials. There are however some differences, and we analyze 
the origin of the difference from a microscopic viewpoint. The band shape of CuA102 is found to 
be even more ideal than that of NaxCo02, and we predict that once a significant amount of holes 
is doped in CuA102, thermoelectric properties (especially the power factor) even better than that 
of NaccCo02 can be expected. 

PACS numbers: 72.15.Jf, 71.20.-b 



I. INTRODUCTION 

Search for good thermoelectric materials serves as an 
intriguing challenge both from the viewpoint of funda- 
mental physics as well as device applications [1]. Ther- 
moelectric materials are often found in semiconductors, 
but the discovery of a large Seebeck coefficient found in 
Naa;Co02[2] which exhibits a metallic nature of the con- 
ductivity opened up a new avenue for the search of ther- 
moelectric materials. Coexistence of the large Seebeck 
coefficient S and the low resistivity p can give rise to a 
large power P = S'^/p, which is important from a applica- 
tion point of view. Soon after the discovery, the origin of 
the large thermopower in Naa;Co02 was studied theoreti- 
cally, where the possible importance of the orbital degen- 
eracy [31,13 or the narrow band width [5[ has been pointed 
out. Later on, one of the present authors along with Arita 
proposed that the peculiar band shape in which the top 
is flat but bends sharply into a dispersive portion plays 
an important role in the coexistence of the large See- 
beck coefficient and the low resistivity in Naa^Co02[6]. 
This band has been referred to as the "pudding-mold" 
type. Recently, various materials with good thermoelec- 
tric properties^uch as CuRhO20-l2|, Li2Rh04[lQ|, |Tlj, 
and Fe As2 [12|, [l3| have been shown to possess this type of 
band shape. The effect of the interplay between electron 
correlation and this kind of band shape has also been 
studied [III, [iMlj . 

Along this line of study, here we focus on CuA102. 
In fact, several previous studies have pointed out a 
strong potential of this material as a good thermoelec- 
tric material pJlflSj. although sufficient carrier doping to 
reduce the resistivity has not been successfully accom- 
plished so far. In ref.[l7, 19], band structure calcula- 
tion has been performed, which shows a flat portion at 
the top of the band, which is reminiscent of the band 
shape of Naa;Co02. In the present paper, we first con- 



struct a model Hamiltonian of CuA102 based on the first 
principles band calculation, and calculate the Seebeck 
coefficient. Then, we compare the model with that of 
Naa^Co02, and discuss the similarities and the differ- 
ences. We also analyze the origin of the difference of 
the band structure between the two materials from a mi- 
croscopic point of view. The pudding-mold type band 
shape of CuA102 is found to be even more ideal than 
that of Naa:Co02, and we predict that once a signifi- 
cant amount of holes is doped in CUAIO2, thermoelectric 
properties (especially the power factor) even better than 
that of Naa^Co02 can be expected. This is due to the 
combination of an extremely ideal band shape (enhances 
the Seebeck coefficient) and a relatively large band width 
(can give large conductivity when holes are doped). 



II. CONSTRUCTION OF THE SINGLE 
ORBITAL MODEL 

We perform first principles band calculation of CUAIO2 
using the Wien2k package |2Q]. We adopt the lattice 
structure parameters given in ref.[19]. Here we take 
RKma^ = 7, 512 /c-points, and adopt the PBE-GGA ex- 
change correlation functional [21]. The calculation result 
is shown as dashed lines in FiglSl which is essent ially 
the same as those obtained in previous studies [17|, |l9|. 
Namely, there is a flat portion around the Brillouin zone 
edge, which bends sharply into a dispersive portion as the 
r point is approached. This is nothing but the pudding- 
mold type band introduced in ref.[6]. 

Since the flat portion of the band is well isolated from 
other portions, we can extract this portion and construct 
a single orbital model, where we exploit the maximally 
localized Wannier orbitals[22, 23]. Here, the extracted 
orbital is obtained by projecting onto the 3(^3^2 _^2 or- 
bital, namely, the obtained Wannier orbital has strong 





FIG. 1. The lattice structure of CuA102(left) and 
NaxCo02 (right) 



III. SEEBECK COEFFICIENT AND POWER 
FACTOR 

A. The Boltzmann's equation approach 

In this section we calculate the Seebeck coefficient 
adopting the single orbital model and using the Boltz- 
mann's equation approach, In this approach, tensors of 
the Seebeck coefficient S and the conductivity a are 
given as. 



S=^K,K^' 



e'Ko 



(1) 



(2) 



where e(< 0) is the elementary charge, T is the temper- 
aturetensors Ki,K2 are given as 
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FIG. 2. (a)The Brillouin zone (b)(left) The first principles 
band calculation result of CUAIO2 (dashed lines) and the sin- 
gle orbital model (solid line) (right) the density of states. 



Cu 3(^322 _^2 character, while the contributions from other 
orbitals due to the hybridization are implicitly taken into 
account in this single Wannier orbital. The result is 
shown in Fig|2]as a solid line superposed to the first prin- 
ciples result. 



Here, e(/c) is the band dispersion, v(}z) = ■|Vfce(/c) is the 
group velocity, r{k) is the quasiparticle lifetime, /(e) is 
the Fermi distribution function, and /i is the Fermi level 
(chemical potential). Here, due to the derivative of the 
Fermi distribution function d/(e)/de, large contributions 
to Kq and Ki come from within ksT from the Fermi 
level /i. In the present study, we approximate r as a 
constant, so that it cancels out in the Seebeck coefficient. 
We simply write axx and Sxx as a and 5, respectively, a 
and thus the power factor aS'^ contains the constant r, 
whose absolute value is not determined. Therefore, we 
only discuss the values of the power factor normalized by 
its maximum value as a function of hole doping rate. 



B. The pudding mold type band 

As discussed in ref. [6] , pudding- mold type band is ad- 
vantageous in obtaining large Seebeck coefficient despite 
low resistivity. For a constant r, eq. (|3|) can roughly be 
approximated as 
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where '^above and I'beiow are group velocities above and 
below the Fermi level (representative values within ksT 
from the Fermi level). Since the Seebeck coefficient is 
proportional to Ki/Kq^ a larger difference between '^above 
and I'beiow gives larger S. Physically, this means that a 
large difference in the group velocities of holes and elec- 
trons results in a large Seebeck coefficient. When the 
Fermi level lies near the band edge, this difference (the 
'^beiow/'^above ratio) cau be large but the absolute val- 
ues of the velocities are small, so that the conductivity 



becomes small. On the other hand, in usual metallic sys- 
tems, in which the Fermi level lies in the middle of the 
bands where '^above and '^beiow have similar values, the 
Seebeck coefficient tends to be small. This is the rea- 
son why a large power factor aS'^ is usually difficult to 
obtain. For the pudding-mold type band, however, the 
large density of states at the top of the band prevents the 
Fermi level from going down rapidly even when a large 
amount of carriers (holes in the present case) is doped, 
and when the Fermi level sits close to the bending point 
of the band, Ki is large because of the small '^above due to 
the flat portion and the large 'i^beiow due to the dispersive 
portion of the band. In this manner, the coexistence of 
large Seebeck coefficient and small resistivity is realized 
for a wide range of hole doping ratio. 



C. Calculation Results 

Figl3] shows the Seebeck coefficient calculated for the 
single orbital model of CUAIO2. FiglSl^a) shows the tem- 
perature dependence for 10% hole doping, and FiglSfb) 
is the doping dependence at T = 300K. These calcu- 
lated values (e.g. 150/iV/K at T = 300K for 10% dop- 
ing) are similar to those values calculated for the single 
band model of NaajCo02 in ref.[6]. It is worth mention- 
ing that in ref. ^] , the band width was reduced to ::^ leV 
(by hand) so as to fit the angle resolved photoemission 
data, while in the present case the band width is about 
4eV, i.e., four times larger. A larger band width implies 
a larger group velocity in the dispersive portion of the 
band, which can give rise to better conductivity. The See- 
beck coefficient being nearly the same despite the much 
larger band width implies that the band shape of CUAIO2 
is even more ideal than that of NaajCo02 from the view- 
point of obtaining large thermopower. Namely, the group 
velocity of the electrons is ideally close to zero (the band 
above the Fermi level is very flat), so that only the holes 
contribute to the thermopower. 

The maximum power factor is reached around 10 per- 
cent hole doping, which is a rather large dopin^rate, and 
is in fact similar to the situation in Na^Co02[6|. This is 
another feature peculiar to the pudding-mold type band, 
where a large amount of doping does not result in a rapid 
reduction of the Fermi level. 



IV. ORIGIN OF THE PUDDING-MOLD TYPE 
BAND 
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FIG. 3. The calculated thermoelectric properties of CUAIO2. 
(a) The Seebeck coefficient against temperature for 10% hole 
doping (b) The Seebeck coefficient against the hole doping 
concentration at T = 300K. (c) The normalized power factor 
against the doping concentration at T = 300K. 



A. The effect of the 2nd and the 3rd nearest 
neighbor hoppings 

Both CUAIO2 and Naa;Co02 exhibit pudding- mold 
type band. Here we discuss its origin from the viewpoint 
of the hopping integrals on the triangular lattice. Here 
we focus on the band structure within the planes, and 
neglect the hopping integral in the z direction. In FigHl 



we show the band structure of the single band model 
of both materials in the Brillouin zone of the two di- 
mensional triangular lattice, in which the hoppings up to 
fifth nearest neighbors are extracted (ti ~ ^5, the orig- 
inal model contains very small hoppings between more 
distant sites). The band structure is normalized by the 
nearest neighbor hopping ti. Note that the band struc- 



TABLE I. The hoppings for the single 
CUAIO2 and Na^Co02 


orbital model of 




ti[eV] ts/ti ts/ti 


tA/tl h/ti 


CuAlOs 
Na^Co02 


-0.45 0.003 -0.069 
0.18 -0.24 -0.15 


-0.030 0.034 
-0.001 0.008 



ture of CUAIO2 is turned up-side down compared to the 
original one because ti < 0. The hoppings normalized 
by ti are given in table [H From this comparison, it can 
be seen that the flat portion of the band occurs around 
the Brillouin zone edge in CUAIO2, while the flat por- 
tion of the band structure of Naa;Co02 is around the T 
point. Note that the path along the Brillouin zone edge 
-K-M- in the two dimensional Brillouin zone corresponds 
to lines like -a-L(or F)- in the three dimensional Brillouin 
zone shown in Figl2fa). 

To see the origin of this difference between the two 
materials, we vary the hopping integrals ^2 and ts "by 
hand". As shown in Fig|5l the second nearest neighbor 
hopping has a dramatic effect of making the band around 
the T point flat while making it around K-M dispersive. 
Consequently, the peak structure of the density of states 
moves largely from £^ < to £^ > as the absolute value 
of ^2 < is increased. The third nearest neighbor hop- 
ping ^3 < on the other hand has the effect of reversing 
the dispersion around K-M, so that the band around K-M 
once becomes nearly perfectly flat (around ts/ti = —0.1) 
as 1^3 1 is increased. 

Now, if we go back to the hopping integrals of the two 
materials given in tabl^ there is a large difference, i.e., 
1^2/^1 1 is almost negligible in CUAIO2 compared to that 
in Naa;Co02. This, along with ts/ti being somewhat 
close to —0.1, makes the band around K-M fairly flat in 
CUAIO2. Strictly speaking, ts/ti ~ —0.07 is still some- 
what away from ts/ti = —0.1, where a nearly perfectly 
fiat band appears. Actually, we find that there is also the 
effect of ts/ti > 0, which makes the band around K-M 
fiat, while enhancing the dispersion around the F point. 
In sharp contrast to CUAIO2, the combined effect of ^2 
and ts in Naa;Co02 makes the band around the F point 
fiat, while enhancing the dispersion around the Brillouin 
zone edge K-M. 



B. The origin of the difference between CUAIO2 
and NaccCo02 

In this section, we discuss the origin of the difference 
in the second neighbor hopping between CUAIO2 and 
NaajCo02. To see this, we now construct models of these 
materials where all of the five 3d orbitals are considered 
explicitly. Here, the difference between the d^z'^-r'^ or- 
bital in the five orbital model and the d^z'^-r'^ orbital 
mentioned in the previous sections should be noted. The 
dsz'^-r'^ Wannier orbital in the previous sections consists 
not only of the d^z'^-r'^ in the five orbital sense, but also 
of the other hybridized orbitals as well. In other words. 




K M K 



FIG. 4. (a)The definition of the hoppings on the triangular 
lattice (b) the two dimensional Brillouin zone, (c)(d) band 
dispersion normalized by ti and the density of states of the 
two dimensional models of CUAIO2 and NaccCo02 which con- 
siders the extracted hoppings ti ^ t^. 
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FIG. 5. Band dispersion of the single orbital model on the 
triangular lattice in which only ti and t2(a) or t3(b) are con- 
sidered. 



the dsz'^_r2 Wannier orbital of the single orbital model 
considered in the previous sections takes into account the 
effect of other orbitals implicitly. Since the band struc- 
ture in the vicinity of the Fermi level is the same between 
the single and the five orbital models, the two models give 
the same transport properties. 

In the five orbital model, the most relevant band orig- 
inates from the d^z'^-r'^ orbital |24| in both CUAIO2 and 
Naa;Co02, but other orbitals can have contribution to 
the band shape due to the hybridization. To see this ef- 
fect, we first hypothetically vary the on-site energy of the 
orbitals other than d^z'^-r^- We find that the d^z/yz oi"- 
bitals have small effect on the band shape in both of the 
materials. On the other hand, we find that varying the 
on-site energy of the d^.y/x'^-y'^ orbitals affects the fiat- 
ness of the band top only in CUAIO2. This means that 
there is a hybridization between d^z2_^2 and dxy/x'^_y2 
orbitals which plays an important role in producing the 
pudding mold type band in CUAIO2. 

To further investigate this point, we consider the ef- 
fects of the intra and interorbital hoppings step by step. 
In FigEl^a) and (c), we show the band dispersion in which 
only the d^z'^-r'^ orbital is extracted from the five 3d 
orbitals, and only the nearest neighbor (a) or the sec- 
ond nearest neighbor hopping (c) is considered. In the 
band dispersion shown in FiglHl^b), we extract three or- 
bitals out of five, namely, d^z'^-r'^ and <ia;y/a?2_^2, and 
consider only the nearest neighbor interorbital hoppings 
between dsz'^_r2 and d^y/x'^-y'^ orbitals. Here, the thick- 
ness of the lines represents the magnitude of the d^z'^_r2 
orbital weight. This result shows that the effect of the 
interorbital d^z'^_r2-dxy/x2_y2 orbital hopping pushes up 
the dsz'2-r'^ band top around the M point. This ef- 
fect is taken into account as a positive contribution to 
12/ ti in the single orbital model which considers the ef- 
fect of d^y/^2_y2 implicitly. Namely, the hopping path ^.j^^ difference in the ^3 



is large and the fiatness of the band does not appear 
around K-M-K. 

The origin of the difference between the two materials 
in the contribution of d^z'^-r'^ -^ d^y/x'^-y'^ ~^ dsz'^-r'^ 
path comes from the difference in the energy level offset 
between d^z-^-r"^ and dj.yij.2_y2 in the five orbital model 
/S.Ed = E{dj.y/j.2_y2) — E{dsz-2-r'^)- As shown in Fig III 
this offset is large in Naa;Co02, while relatively small in 
CUAIO2. A smaller energy difference gives larger con- 
tribution of the dsz'2-r'^ -^ djy/j2_y2 -^ d^z'^_r2 path. 
The difference in AE^ can be understood from the lat- 
tice structure. As shown in the bottom of Figd the 
oxygen atoms in CUAIO2 are located at positions toward 
which the d^z'^_r-2 orbitals are elongated, while that is 
not the case for Naa^Co02. Therefore, the crystal field of 
the ligand atoms pushes up the d^z'^-r'^ level, locating it 
just above the djy/j2_y2 level. 



V. CONCLUSION 

To conclude, CUAIO2 is a very good candidate for ther- 
moelectric material with large Seebeck coefficient coex- 
isting with large conductivity. The origin of this is the 
pudding mold band similar to, but different from, the 
one in Naa;Co02. In CUAIO2, the negligibly small sec- 
ond nearest neighbor hopping does not affect the fiat por- 
tion of the band around the Brillouin zone edge already 
present in the nearest neighbor hopping model on the tri- 
angular lattice. Moreover, the additional presence of the 
third (and the fifth) hopping integrals makes the band 
even more flat around the Brillouin zone edge. This is 
in sharp contrast with the case of Naa;Co02, where the 
large second nearest neighbor hopping moves the flat por- 
tion of the band to the F point area. The difference of 
the band shape between the two materials comes from 



dsz'^-r'^ -^ dxy/x'^-y'^ ~^ dsz'^-r'^ givcs a positivc con- 
tribution to the effective ^2/^1 between second neighbor 
(i3^2_^2 orbitals. On the other hand, there is a negative 
contribution to ^2/^1 of the single orbital model coming 
from the direct second nearest neighbor d^z'^-r'^-dsz'^-r'^ 
hopping in the flve orbital model, which pushes down 
the band around the M point, as shown in FigEl^c). In 
CUAIO2, the positive and negative contributions almost 
cancel with each other, resulting in a very small ^2/^1 in 
the single orbital model. 

A similar flve orbital analysis for Naa;Co02, shown in 
Figsl6l^d)-(f), reveals that the situation is different. From 
FigEl^e), it can be seen that the dispersion originating 
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hybridization is very small. 



namely, the contribution of the d^z'^_j.2 -^ djy/j2_y2 -^ 
d^z2_r2 path to ^2 in the single orbital model is negligi- 
ble. On the other hand, there is a direct (i322_^2-(i322_^2 
second neighbor hopping, which gives a positive contri- 
bution to 12/ ti in the single orbital model, pushing up 
the band around the M point. Hence, in Naa;Co02, ^2/^1 



path in the flve orbital sense ; a fairly close energy lev- 
els between d^z'^_r2 and djy/j2_y2 gives rise to a large 
contribution from this path, which cancels out the di- 
rect dsz'2-r'^-dsz'2-r'2 contribution to the second neighbor 
hopping. The difference in the energy difference between 
d^z'^_r2 and djy/j2_y2 levels comes from the position of 
the oxygen atoms. A large Seebeck coefficient compara- 
ble to the model of Naa^Co02, despite a much larger band 
width, implies that the band shape is extremely ideal, 
and a very good thermoelectric properties, especially a 
large power factor, are expected once large amount of 
hole doping is accomplished. 
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FIG. 6. Band dispersion of models constructed by extracting some of the orbitals and hoppings from the 3d five orbital models 
of CuA102((a)-(c)) and Naa;Co02((d)-(e)) The thickness of the lines represent the strength of the d^z'^-r'^ orbital character. 



(a)(d) Only the d^^2_r2 orbital extracted, only the nearest neighbor hopping considered. (b)(e) 



and i 



extracted, only the nearest neighbor interorbital hopping considered. (c)(f) Only the d^^2_ 
nearest neighbor hopping considered. 



^xy I x^—y 



orbitals 



orbital extracted, only the second 



[1] For a general review on the theoretical aspects as well 
as experimental observations of thermopower, see, G.D. 
Mahan Good Thermoelectrics, Solid State Physics 61, 81 
(1997). 

[2] I. Terasaki, Y. Sasago and K. Uchinokura, Phys. Rev. B 
56 R12685 (1997). 

[3] W. Koshibae, K. Tsutsui and S. Maekawa, Phys. Rev. B 
62 6869 (2000). 

[4] W. Koshibae and S. Maekawa, Phys. Rev. Lett. 87, 
236603 (2001). 

[5] D.J. Singh, Phys. Rev. B 61, 13397 (2000). 

[6] K.Kuroki and R.Arita, J.Phys. Soc.Jpn. 76 083707 
(2007). 

[7] H. Kuriyama, M. Nohara, T. Sasagawa, K. Takubo, T. 
Mizokawa, K. Kimura, and H. Takagi, in Proc. 25th In- 
ternational Conference on Thermoelectrics (IEEE, Pis- 
cataway, 2006). 

[8] S.Shibasaki, W.Kobayashi, and I. Terasaki: Phys. Rev. B 
74, 235110 (2006). 

[9] H. Usui, R. Arita, and K. Kuroki, Journal of Physics. 
Condensed Matter : an Institute of Physics Journal 21, 



064223 (2009). 
[10] Y. Okamoto, S. Niitaka, M. Uchida, T. Waki, M. Taki- 

gawa, Y. Nakatsu, A. Sekiyama, S. Suga, R. Arita, and 

H. Takagi, Phys. Rev. Lett. 101, 086404 (2008). 
[11] R. Arita, K. Kuroki, K. Held, A. V. Lukoyanov, S. Sko- 

rnyakov, and V. I. Anisimov, Phys. Rev. B 78, 115121 

(2008). 
[12] P. Sun, N. Oeschler, S. Johnsen, B.B. Iversen, and F. 

Steglich, Appl. Phys. Express 2 091102 (2009). 
[13] H. Usui, K. Kuroki, S. Nakano, K. Kudo, and M. Nohara, 

l arXiv:12TL 7176 
[14] P. Wissgott, A. Toschi, H. Usui, K. Kuroki, and K. Held, 

Phys. Rev. B 82, 201106 (2010). 
[15] P. Wissgott, A. Toschi, G. Sangiovanni, and K. Held 

Phys. Rev. B 84, 085129 (2011). 
[16] H. Usui, PhD Thesis, 2012, Univ. of Electro- 
Communications, unpublished. 
[17] H. Katayama-Yoshida, T. Koyanagi, H. Funashima, H. 

Harima, A. Yanase, Solid State Comm.126, 135 (2003). 
[18] P. Poopanya, A. Yangthaisong, C. Rattanapun, A. 

Wichainchai, J. of Electronic Mat., 40 987 (2011). 




dyzhx 

0.2eV 


7^ 

ja 


■ (iiz2-,2 ^^ 


|A£'d|=0.2eV 


Ji^ 



Uxylx^-l 



[19] A. Nakanishi, H. Katayama-Yoshida, Solid State 
Comm.152, 24 (2012). 

P. Blaha, K. Schwarz, G.K.H. Madsen, D. Kvasnicka, 
and J. Luitz, Wien2k: An Augmented Plane Wave + 
Local Orbitals Program for Calculating Crystal Properties 
(Vienna University of Technology, Wien, 2001). 
J. P. Perdew , K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

[22] N. Marzari and D. Vanderbilt, Phys. Rev. B 56, 12847 
(1997); I. Souza, N. Marzari, and D. Vanderbilt, Phys. 
Rev. B 65, 035109 (2001). The Wannier functions are 
generated by the code developed by A. A. Mostofi, J. 
R. Yates, N. Marzari^ L Souza, and D. Vanderbilt, 



[20] 



[21] 



NaxCo02 



CuA102 



(http://www.wannier.org/ ). 

[23] J. Kunes, R. Arita, P. Wissgott, A. Toschi, H. Ikeda, and 
K. Held, Comp. Phys. Commun. 181 1888 (2010). 

[24] For NaccCo02, the following notation is more commonly 
used [5]. The 3d orbitals consist of the low- lying t2g and 
the high-lying Cg orbitals, and the t2g orbitals are split 
into aig and e'g orbitals. In the present notation, the 
d^z2_r2 orbital corresponds to the commonly used aig or- 
bital. We will call this the ''d^z2_r2 orbital" in the sense 
that it is the maximally localized Wannier orbital ob- 
tained by projecting onto the d^z'^_r2 orbital. The bond- 
ing and antibonding states between the present d^y/x2_y2 
and dxz/yz orbitals correspond to the Cg and Cg orbitals, 
respectively. 
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